function diff = solve_b_minimal(D, v2, param, glob, options)
    
    param.D = D;
        
    c = [v2; v2];
    
    % Solve the incumbent's value function
    v     = solve_valfunc(c,glob.s,param,glob,options);
    v1    = v.v1;
    
    % Find the expected value of an entrant
    diff = v1(glob.s(:,1) == 0, :)'*glob.a_stat_E' - param.ce;
        
    
end